15006A Xi V 2HE 193 Avil 2015 


e 
. 


arXiv 


PREPARED FOR SUBMISSION TO JHEP 


Interpreting the galactic center gamma-ray excess in 
the NMSSM 


Junjie Cao,®? Liangliang Shang,“ Peiwen Wu,° Jin Min Yang,° Yang Zhang® 
Department of Physics, Henan Normal University, Xinziang 453007, China 

è Department of Applied Physics, Xi’an Jiaotong University, Xi’an 710049, China 

“State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Academia Sinica, 
Beijing 100190, China 
E-mail: junjiec@itp.ac.cn, shlwel11988@gmail.com, pwwu@itp.ac.cn, 
jmyang@itp.ac.cn, zhangyang@itp.ac.cn 


ABSTRACT: In the Next-to-Minimal Supersymmetric Standard Model (NMSSM), all 
singlet-dominated particles including one neutralino, one CP-odd Higgs boson and one 
CP-even Higgs boson can be simultaneously lighter than about 100 GeV. Consequently, 
dark matter (DM) in the NMSSM can annihilate into multiple final states to explain the 
galactic center gamma-ray excess (GCE). In this work we take into account the foreground 
and background uncertainties for the GCE and investigate these explanations. We carry 
out a sophisticated scan over the NMSSM parameter space by considering various experi- 
mental constraints such as the Higgs data, B-physics observables, DM relic density, LUX 
experiment and the dSphs constraints. Then for each surviving parameter point we per- 
form a fit to the GCE spectrum by using the correlation matrix that incorporates both 
the statistical and systematic uncertainties of the measured excess. After examining the 
properties of the obtained GCE solutions, we conclude that the GCE can be well explained 
by the pure annihilations ¥?v? > bb and XX > Ai Hi; with A, being the lighter singlet- 
dominated CP-odd Higgs boson and H; denoting the singlet-dominated CP-even Higgs 
boson or SM-like Higgs boson, and it can also be explained by the mixed annihilation 
xix? —> WtW-, Ai. Among these annihilation channels, 9%? —> AiH; can provide 
the best interpretation with the corresponding p-value reaching 0.55. We also discuss to 
what extent the future DM direct detection experiments can explore the GCE solutions 
and conclude that the XENON-1T experiment is very promising in testing nearly all the 
solutions. 
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1 Introduction 


The compelling evidences for the existence of Dark Matter (DM) from various cosmological 
and astrophysical observations have provided us a good portal in the search for new physics 
beyond the Standard Model (SM). One possible method to explore DM in the present Uni- 
verse is the indirect detection, which looks for the particles produced when DM annihilates 
in the DM halo. These particles include photons, antiparticles and neutrinos, and among 
them gamma rays have often been defined as the golden channel for DM indirect detec- 
tion since the signal can be traced back to the source. The Large Area Telescope (LAT) 
onboard the Fermi Gamma-ray Space Telescope, due to its unprecedented angular and 
energy resolutions, has produced the most detailed maps of the gamma ray sky for a wide 
range of energies. Intriguingly, as was reported by several independent groups [1-9] and 
also by Fermi Collaboration itself [10], the Fermi-LAT data have revealed the presence of 
an extended excess of gamma rays over the modeled foreground and background emissions 
towards the Galactic Center (GC). Although several astrophysical mechanisms, such as the 
thousands of unresolved millisecond pulsars [11-13] and the interactions between comic rays 
(CR) and interstellar gases [14-18], have been proposed to interpret this Galactic Center 
Excess (GCE), they usually fail to generate the morphology and energy spectrum of the 
GCE simultaneously!. So in this work, we instead consider another possibility that the 
GCE is produced by the annihilation of DM. Although this interpretation has been con- 
strained by the measurements of CR such as the Fermi-LAT detection of the gamma-rays 
from dwarf spheroidal galaxies (dSphs) [20-23], the non-observation of spectral features 
in the AMS-02 measurements of CR positron [24-27], and PAMELA observations of the 
CR anti-protons [28-34], it still remains a most attractive one not only because the excess 
emission shows spectral and morphological properties consistent with a telltale sign from 
DM annihilation, but also because in such an interpretation, the annihilation cross section 
required to explain the GCE is of the right size to account for the DM density from thermal 
freeze-out. 

So far there have been a large number of attempts to explain the GCE by DM anni- 
hilation in various new physics models [35-102]. In the early analyses of the annihilations, 
great efforts were focused on the channels ~y — bb with mg ~ 35GeV and XX > 77 with 
mg ~ 10GeV since they can reproduce well the GCE spectrum obtained at that time. 
Recently a critical reassessment of the DM interpretation was made by examining in a 
comprehensive way the foreground and background uncertainties [9]. It was found that 
taking the estimated uncertainty in the high-energy tail of the spectrum into account, a 
much larger number of DM annihilations are able to fit well the y-ray data than previously 
noted [73, 79]. Explicitly speaking, as far as the annihilation YY — bb is concerned, now 
the mass of DM is extended to a broader range from 30GeV to 70GeV in explaining the 
GCE [73, 79]. Other annihilation channels such as DM annihilation into light quark pairs 
and even gluon pair are also able to provide a good fit to the GCE [79]. More strikingly, 
this new analysis opens up a very good solution usually neglected before, namely DM an- 
nihilation into a pair of light non-standard Higgs bosons [49, 50]. This important progress 


‘An exception may be the mechanisms recently proposed in [19]. 


motivates us to renew the solutions to the GCE in supersymmetric theories, which usually 
predict the lightest neutralino Y? as a natural DM candidate. 

As the most economical realization of supersymmetry, the Minimal Supersymmetric 
Standard Model (MSSM) is unsatisfactory in explaining the GCE due to the following four 
reasons [66, 93]. First, the relic density of DM has required its mass to be larger than 
about 40GeV [103]. In this case, the annihilations ¥9X¥? > 77, qq with q denoting a light 
quark can not provide a good fit any more. Second, except for excessive fine-tuning cases 
the LHC experiments have pushed the lower mass bounds for the CP-odd Higgs boson 
and the bottom squarks up to several hundred GeV. As a result, the cross section of 
DM annihilation into bb in present day is too small to significantly contribute to the GCE 
(66, 93]. Third, due to the small velocity of DM in our galaxy, the annihilation rate for 
DM into SM-like Higgs pairs is p-wave suppressed. Consequently this channel is not large 
enough to generate the GCE. Finally, as for the annihilations 9%} > WW, ZZ, their fits 
to the GCE spectrum indicate that regardless of their annihilation rates the corresponding 
p-values are always less than 0.04 [73, 79]. This means that the annihilations can not 
generate the proper spectrum shape for the GCE. We note that for a given parameter 
point of the MSSM, DM usually annihilates into multiple final states. In this case, the 
situation can not be improved greatly because, due to the particle spectrum of the MSSM 
allowed by the current experiments, either the total cross section falls short for the GCE, 
or the dominant annihilation channel can not reproduce the GCE spectrum well [92]. 

Given the problems of the MSSM, we consider to interpret the GCE in the Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) with a Z3 symmetry, which is the 
simplest gauge singlet Higgs extension of the MSSM [104]. Distinguished from the MSSM, 
the NMSSM predicts three singlet-dominated particles: one neutralino, one CP-even and 
one CP-odd Higgs bosons. These particles are rather special in that all of them can be 
simultaneously lighter than about 100GeV, and that the couplings for the interactions 
among themselves are determined by the parameter Kk, which alone is able to predict the 
right rates for some annihilation channels to explain the GCE (see the following discussion). 
These features make the NMSSM with a singlet-dominated DM well suit to account for 
the GCE because, as we will show below, some golden channels for the GCE need light 
particles to act as the DM, the mediator and/or the annihilation final state. 

We note that the interpretations of the GCE in the NMSSM have been intensively 
discussed in [57, 58, 66, 68, 69, 93]. However, in [57, 58, 66, 68, 69] the authors did not 
consider the systematic uncertainties mentioned above. As a result, the model parameter 
space they considered is much narrower than that of this work and the obtained conclusions 
were incomplete. While for [93], although the authors have taken the uncertainties into 
account, they considered the parameter space characterized by a large A which is different 
from our discussion. 

The aim of this work is to explore any possible solution to the GCE in the Z3 NMSSM. 
For this end, we perform a sophisticated scan over the model parameters by considering 
various experimental constraints such as the DM relic density, the Higgs data as well as 
the observation of dwarf galaxies. We use the correlation matrix presented in [9] to include 
the systematic uncertainties on the GCE spectrum and only keep the parameter points 


that can reproduce well the spectrum. In our study we mainly consider a singlino-like DM 
which is believed to interpret the GCE without excessive fine tuning. As we will show 
below, the annihilation ¥?v? —> H;Aı with H;A, denoting a scalar-pseudoscalar Higgs pair 
may provide the best fit to GCE, and the canonical annihilation X9? — bb still remains 
a satisfactory solution except that mgo is now allowed to vary within a broader range. 
Moreover, it is interesting to see that the mixed annihilation into W*W~— and H;A; final 
states is also able to generate a spectrum consistent with the GCE. These conclusions are 
quite different from previous studies in the NMSSM. 

This paper is organized as follows. In Section II, we introduce some of the characteristic 
features of NMSSM, the basic knowledge about the GCE and our strategy for the parameter 
scan. In Section III, we discuss in detail the interpretations of the GCE when Hə is the 
SM-like Higgs boson, and in Section IV, we carry out a similar study but for the case that 
H; acts as the SM-like Higgs boson. We draw our conclusion in Section V and provide 
more information of the NMSSM couplings in the Appendix. 


2 Fitting the GCE in the NMSSM 


2.1 Theoretical setup for the GCE in the NMSSM 


We start our analysis by recapitulating the basics of the NMSSM. As one of the most eco- 
nomical extensions of the MSSM, the NMSSM introduces one gauge singlet Higgs superfield 
in its matter content, and since one purpose of the extension is to solve the y-problem of 
the MSSM, a Z3 symmetry is usually adopted in the construction of the superpotential to 
avoid the appearance of parameters with mass dimension. As a result, the superpotential 
of the NMSSM and the soft breaking terms in Higgs sector are given by [104] 
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where Wp is the superpotential of the MSSM without the p-term, H, Hy and $ are Higgs 
superfields with H,,, Hq and S acting as their scalar components respectively, the dimen- 
sionless coefficients À and « parameterize the strengthes of the Higgs self couplings, and 
Mu, Ma, Mg, A, and A, are soft-breaking parameters. In practice, after the electroweak 
symmetry breaking the soft-breaking squared masses 72, m3 and mn? are traded for mz, 
tan 8 = v,/vg and u = Avs as theoretical inputs. 

Due to the presence of the superfield $, the NMSSM contains a singlino field which is 
the fermion component of S$, and one more complex Higgs field S compared to the MSSM. 
As a result, the neutralino mass eigenstates Y? (with 7 ranging from 1 to 5) are the mixtures 
of bino, wino, higgsinos and singlino, and the CP-even (odd) Higgs mass eigenstates H; 
with i = 1,2,3 (A; with i = 1,2) are mixtures of the real (imaginary) parts of H,, Ha 
and S. Throughout this paper, we assume the mass order mgo < Mgo < +++ < mgo for 
neutralinos, and my, < MH, < MHz, MA, < MA, for Higgs bosons. 

There are three distinguished features in the NMSSM. One is that DM in the NMSSM 


may be either singlino-dominated or bino-dominated. As expected, the properties of a 


singlino-dominated DM are quite different from those of a bino-dominated DM, which 
makes the DM physics in the NMSSM much richer than that in the MSSM [105]. Another 
feature is that, in the presence of a singlino-dominated DM with mass below 100GeV, the 
singlet-dominated CP-even and CP-odd Higgs bosons can be simultaneously lighter than 
about 100GeV [105, 106], and the strengthes for the interactions among these particles are 
determined by the parameter « which may be as large as 0.1. This feature, as we will show 
below, makes the NMSSM with a singlino-dominated DM well suit to explain the GCE. 
In the appendix, we list the properties of these particles used in our analysis. The other 
feature is that either Hı or Hə in the NMSSM can act as the SM-like Higgs boson [107] 
and generally speaking, Hz as the SM-like boson is more attractive from phenomenological 
point of view and also from naturalness argument. 

In the DM explanation of the GCE, the observed y-ray originates mainly from the 
cascade decays of the annihilation final states. In the NMSSM, the possible annihilation 
final states include ff, VV, H;H;, A;A; and H;A; [108], where f (V) denotes any of the 
fermions (vector bosons) in the SM, and H; (A;) denotes a CP-even (CP-odd) Higgs boson. 
In this work, we are particularly interested in the annihilations %9} > bb, WtW, Hj A. 
These annihilations proceed through the s-channel mediator of a Z boson or a Higgs boson 
with an appropriate CP quantum number, and also proceed through the t/u-channel ex- 
change of a sbottom, a chargino and a neutralino respectively. The complete expressions 
of the annihilation cross sections are rather complicated, but in non-realistic limit, i.e. the 
velocity of DM approaching zero, some contributions become unimportant. In this case, 
the velocity weighted annihilation cross section can be approximated by [108] 
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where C'xyz denotes the coupling of the interaction involving the particles X, Y and Z, 
T4; is the width of the CP-odd state A;, w = mo /miy, ki = me. /miy, fit (fir) is 
1 


i 


the coupling coefficient for KX LW? (XIX RWT) interaction, and ô = (2mgo — (my, + 
ma,))/2mzo. In getting Eq.(2.5), we note that a good fit to the GCE requires that the 
HA, final state is produced close to threshold, i.e. 6 ~ 0, so we can expand (o4,4,V)0 
in terms of ô. Then the first two terms on the right hand of Eq.(2.5) come from the left 
diagram of Fig.1, and the last term comes from the right diagram of Fig.1. 

The flux per unit solid angle at some photon energy Ey, which is observed by Fermi- 


Figure 1. Feynman diagrams contributing to the annihilation 7x? > H:Aı with X$ (j from 1 to 
5) denoting any of the five neutralinos. A u-channel diagram in associated with the t-channel one 
is assumed. 


LAT, is then given by 
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where dNj.,-/dE is the photon spectrum generated by the annihilation xoxo > XY, ppm 
is the DM profile and the integral over p},, is along the light-of-sight (LOS) at an angle 0 
towards GC. In the DM interpretation of the GCE, a generalized Navarro, Frenk & White 
(NFW) DM profile is usually adopted, and its expression is given by [109, 110] 


plr) = po Op HH T (2.7) 


with slope parameter y = 1.26, scale radius Rs = 20 kpc and the local DM density 
po = 0.4GeV/cm? at the radial distance of the sun from the galactic centre ro. Here 
the coordinate r is centered on the galactic centre and can be expressed as r?(s,0) = 
rå + s — 2ros cos with s and @ being the LOS distance and the aperture angle between 
the axis connecting the earth with the galactic centre and the LOS respectively. 

In our study, we use the package micrOMEGAs-3.6.9.2 [111] to calculate the DM relic 
density and with the help of PYTHIA [112] to generate the flux in Eq.(2.6). Note that 
in any explicit model, DM usually annihilates into multiple final states. In this case, the 
different fluxes are summed over. 


2.2 Parameter scan strategy for GCE solution 


We simplify our scan over the NMSSM parameter space by fixing the parameters that are 
not closely related to the DM studies. The soft SUSY breaking parameters in the squark 
sector are all fixed to be 2 TeV except that we vary those for the third generation to 
generate a CP-even Higgs near 125 GeV. We assume A; = A, and My, = Mp, to reduce 
the number of free parameters. Similarly, all of the soft SUSY breaking parameters in 
the slepton sector are fixed to be 300 GeV to explain the discrepancy of the measured 
value for muon anomalous magnetic moment from its SM prediction. As for the gaugino 
sector we abandon the Grand Unified Theory relation and fix the wino mass and gluino 
mass at 1 TeV and 2 TeV respectively. Consequently, the remained free parameters include 
tan p, u, À, K, A), Ag in the Higgs sector, Mg,, My, and A; for third generation quarks 


and the bino mass Mj, which are all defined at the scale of 2 TeV in the scan. We use 
NMSSM Tools-4.3.0 [113] to scan intensively the following NMSSM parameter region: 
1<tanB < 40, 0<A<0.7, 0< |K| < 0.7, |Mi| < 600 GeV, 
0< A, <5 TeV, |A,| <2 TeV, |A| < 5 TeV, 
100 GeV < u < 600 GeV, 200 GeV < mQ, Mu; < 5 TeV. (2.8) 


Table 1. Favored parameter region of the NMSSM to explain the GCE, which are classified by the 
These annihilations are called Solution I, H, III, IV and 
V respectively in the following discussion. All input parameters are defined at 2TeV and quantities 


dominant final state in DM annihilations. 


with mass (annihilation cross section) dimension are in unit of GeV (10 


-26 em? /s). 


Hə is SM-like Hı is SM-like 

bb AYA, wtw- ALA, H Aı 
tanB (10, 30) 8, 40) 8, 40) 15, 20) 12,20) 
Mı (—600, —60) —600, —80) —600, — 130) —130, —90) —200, —100) 
AÀ (0.2, 0.7) 0.2, 0.6) 0.2, 0.4) 0.6, 0.7) 0.4, 0.7) 
K (0.02, 0.12) 0.07, 0.15) 0.09, 0.14) 0.10, 0.14 0.11, 0.16) 
H (160, 300) 110, 210) 110, 160) 220, 270) 210, 270) 
Ay (2400, 5000) 830, 5000) 970, 5000) 3900, 5000) 2900, 5000) 
Ax (—210, —70) —60, 22) —70, 10) —65, —16 —66, 5) 
Atb (—4300, 3900) —4600, 4700) —4700, 3900) —2200, 2000) —3400, 4000) 
Mas; (300, 5000) 350, 5000) 500, 5000) 1200, 4600) 700, 4800) 
Mvz,p3 | (250, 5000) 270, 5000) 400, 5000) 250, 5000 1400, 5000) 
MH (15, 102) 61,119) 83, 110) 124, 127.4) 124.5, 127.4) 
MH» (122.8, 127.8) 122.7, 128) 123, 128) 125.6, 142) 125.7, 146) 
(ov)o (0.17,1.9) 0.29,1.8) 0.44,1.6) 0.34,1.2) 0.38,1.5) 
mgo (31,70) 62, 114) 84, 102) 71, 87) 80, 127) 
M0 (67, 298) 83, 233) 114, 165) 86, 128) (100, 192) 
m (166, 297) 117,214) 117, 158) 225, 270) (218, 266) 
MA, (58, 133) 9, 109) 10, 105) 9, 38) (16, 95) 
My (3510, 4666) 1477, 3447) 2818, 2968) 3883, 4941) (2928, 4740) 
Xe vin (23, 35) 21, 35) 24, 35) 21, 35) (21, 35) 


The process to retain the parameter points include the following steps: 


e We require the DM to be singlino-dominated and satisfy mgo < 150 GeV, and impose 
all the experimental constraints encoded in NMSSMTools-4.3.0 [113] which include 
the relic abundance at 3ø level (0.107 < Qh? < 0.131), LUX exclusion bound at 90% 
C.L., 
moment at 20 level. We also consider various electroweak precision data calculated 
n [114]. 


various B-physics measurements as well as the discrepancy of muon magnetic 


e We consider the constraints on the Higgs sector with the package HiggsBounds-4.1.2 
[115] which contains the data from LEP, Tevetron and LHC. For the SM-like Higgs 
boson, we further perform a fit to the data with the package HiggsSignal [116] and 
keep the 2o samples. 


e We use micrOMEGAs-3.6.9.2 [111] to calculate the DM annihilation cross section at 
present day, and then impose the constraints from dSphs by the data in [23] for the 
bb annihilation channel and with the method introduced in [93] for the H;Aı final 
states. 


We also use micrOMEGAs-3.6.9.2 [111] to generate the y-ray spectrum. Considering 
the astrophysical uncertainties which may come from the errors in our setting on the 
local DM density po, the scale radius R, and the inner slope parameter y in Eq.(7), 
for each parameter point we allow an uncertainty factor A in the range of (0.17, 5.3) 
for the annihilation cross section, or equivalently for the height of the gamma-ray 
spectrum in Eq.(2.6) [79]. Then for the A—tuned y-ray spectrum, we perform a 
fit to the residual GCE spectrum obtained in [9] by using the publicly available 
covariance matrix, which include both the statistical and systematic uncertainties of 
the measured flux. The corresponding en function is calculated by [9, 79]: 


dN dN dN dN 

2 -1 

= ee Pe Wd or 2.9 

where X;; is the covariance matrix, dN /dE; is the measured flux in the i-th energy bin, 
and dN /dE; is the flux predicted by the NMSSM, which depends on the parameter 
point and also on the factor A. 


We define the GCE x? as the minimum value of x3,(A) among different choices of 
A, Xécpg = min(x2,(A)), and keep the parameter points that satisfy yacg < 35.2. 
These points are assumed to have the capability to explain the GCE at 95% confidence 
level for 23 degree of freedom [9]. 


The parameter ranges of the GCE solutions are listed in Table.1, which are classified 
by the dominant final state in DM annihilations (see the following discussion). For the first 
three types of the DM annihilations Hə acts as the SM-like Higgs boson, while for the last 
two types Hı corresponds to the SM-like Higgs boson. One distinguished feature that Table 
1 exhibits is that all the singlet dominated particles in the GCE solutions, including DM, 
the singlet-dominated CP-even and CP-odd Higgs bosons, are lighter than about 150GeV. 
This feature, as we will emphasized below, makes the NMSSM well suit for explaining the 
GCE. 


3 GCE solutions with H being the SM-like Higgs boson 


In this section, we exhibit the features of the GCE solutions for the case that DM is 
singlino-dominated and Hə acts as the SM-like Higgs boson. All the solutions considered 
in this work survive the constraints listed in last section and meanwhile can explain the 
GCE at 95% C.L.. 

In Fig.1 we project the solutions on (av)9 — Mo plane (upper panel) and xĝcp — mgo 
plane (lower panel). Solutions marked by red square, blue triangle and black asterisk 


correspond to the cases that DM annihilates with the largest branching ratio into bb, Hy Aı 


Table 2. Detailed information of the benchmark points used in our discussion. Quantities with 
mass, annihilation and scattering cross section dimension are in unit of GeV, cm?/s and pb respec- 


tively. 

Point tanB A K H Ay Ar ADs,Us3 Mı MQ; Mvu;,Ds 
Pl 16 0.36 0.04 241 3891 -136 420 -472 4127 4445 
P2 12 0.46 0.12 179 2036 -6 -2354 -209 2197 3673 
P3 13 0.27 0.11 130 1899 -5 -524 -170 4098 4384 
P4 18 0.69 0.12 243 4518 -43 -320 -103 1436 4308 
P5 17 0.66 0.13 226 3923 -17 1138 -97 4540 1286 
P6 18 0.66 0.15 217 4048 -24 2050 -103 4170 1452 
P7 15 0.50 0.13 255 4085 -35 2621 -131 2935 4468 

Point | ma, | MH, | Myt | Ma, myo | mgg Migs Br (n2-+x9x2) Brini+a,A,) | Br(n2+a, Aj) 
Pl 40 125 3960 99 50 256 248 0.54% 0 0 
P2 99 125 2065 66 87 178 183 0 0 0 
P3 99 126 1823 48 92 126 134 0 88.51% 7.06% 
P4 126 133 4452 20 78 102 249 0 6.50% 95.32% 
P5 125 126 3883 27 81 96 231 0 4.86% 94.01% 
P6 126 129 4022 33 85 101 222 0 4.55% 95.74% 
P7 125 145 4068 69 121 127 262 0 0) 94.11% 

Point | yéce | p-value | (ov)o | (ovire as” ae? Ry Riya, | Rytw- | Rapa, 


Pl 23.3 0.44 6.1E-27 | 2.8E-26 | 2.3E-15 | 1.5E-04 | 90.4% 0.0% 0.0% 0.0% 


P2 21.6 0.54 1.4E-26 | 2.6E-26 | 8.2E-10 | 1.3E-03 | 0.0% 96.8% 3.1% 0.0% 


P3 24.8 0.36 1.3E-26 | 2.7E-26 | 1.9E-10 | 8.5E-04 | 0.1% 46.1% 47.2% 0.0% 


P4 22.4 0.50 8.6E-27 | 2.8E-26 | 5.5E-10 | 1.6E-03 | 0.1% 94.3% 0.1% 5.4% 


P5 21.4 0.55 9.9E-27 | 3.1E-26 | 4.4E-10 | 1.7E-03 | 0.1% 68.3% 0.2% 31.4% 


P6 21.6 0.54 8.3E-27 | 3.2E-26 | 9.8E-10 | 2.0E-03 | 0.1% 42.1% 2.1% 55.5% 


P7 23.7 0.42 7.9E-27 | 3.3E-26 | 1.1E-09 | 5.5E-04 | 0.1% 2.2% 7.2% 84.0% 


and WtW- final states respectively, which hereafter are collectively called Solution I, 
Solution II and Solution III correspondingly. Then the upper panel of Figure.2 indicates 
that, for the ranges 30GeV < mgo < 40GeV, 50GeV < mgo < 62GeV and 63GeV < 
mgo < 70GeV, Solution I is viable, while for 683GeV < mgo < 115GeV and 83GeV < 


mo < 100GeV, Solution II and Solution III can account for the GCE respectively. For any 
of the solutions, the (cv)o is larger than 1.7 x 10~?"cm3/s, and its lower bound increases 
monotonically as ¥? becomes heavier. The reason for the latter behavior is that, for a 
heavier DM, its number density is smaller. So to obtain the same photon flux for the GCE, 
a larger cross section is needed. 

The lower panel of Fig.2 indicates that the best interpretation in Solution I comes 
from mgo ~ 50GeV with oe p œ 23 and a p-value of 0.44. This conclusion coincides 
with that of [79], which was obtained in a model independent way and for a pure bb 
annihilation channel. For Solutions II and Solutions III, the best interpretations locate at 
mgo ~ 87GeV with XŻCE ~ 21.6 and mgo ~ 92GeV with Moan ~ 24.7 respectively, and 
the corresponding p-values are 0.54 and 0.36. These two solutions, within our knowledge, 


were rarely discussed in previous literatures about the NMSSM. Moreover, we checked 
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Figure 2. The GCE solutions of the NMSSM for a singlino-dominated DM with H> acting as the 
SM-like Higgs boson, which are projected on the (av)o — mgo plane (upper panel) and y2aopn — 
mgo plane (lower panel). Solutions marked by the red square, the blue triangle and the black 
asterisk correspond to the case that DM annihilates in present day mainly by the channels ¥?x? > 
bb, Hı A1, WtW- respectively, which are collectively called Solution I, II and III correspondingly. 


that, in the case of myo ~ mz /2 (myo ~ 62GeV), DM annihilated in early universe mainly 
through a nearly on-shell Z boson (SM-like Higgs boson). Since nowadays this dominant 
annihilation channel is helicity (p-wave) suppressed, (avo can not reach the size required 
for the GCE. 

In Table 2, we present detailed information of three benchmark points P1, P2 and 
P3 for Solution I, II and III respectively. This table indicates that the sole annihilation 
channel 9%? — bb or X9X? —> Hi A can be responsible for the GCE; while for the channel 
Xx? — WtW-, it must mix sizeably with the channel ¥?¥? + HA; to account for the 
GCE. We will return to this issue later. 

In our calculation, we found that the condition on the GCE y? can reduce the number 
of the parameter points that survive the constraints by more than 90%. This implies 
that the GCE has non-trivial requirements on the parameters of the NMSSM, especially it 
suggests that some of the independent parameters may be correlated. Motivated by this 
thought, we study the correlations among the parameters A, K, H, mgo and ma, which 
are important parameters in the interpretation of the GCE and show the corresponding 
results in Fig.3. In the following, we concentrate separately on each kind of the solutions 
and investigate its features. Such a study is helpful to understand the correlations in Fig.3 
and also the properties of the benchmark points listed in Table.2. 
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Figure 3. Same as Fig.2, but showing the correlations of different parameters. 


3.1 Solution I - the bb annihilation channel 


Among the solutions to the GCE, Solution I is the most intensively studied one. After 
considering the systematic uncertainties, one important improvement of Solution I over its 
previous version is that DM mass is now allowed in the range from 30GeV to 70GeV, which 
is much wider than before. 

The key features of Solution I are as follows: 


e The lighter CP-odd Higgs boson is correlated with DM by ma, ~ 2mxo. This corre- 
lation is shown in the upper right panel of Fig.3 which means that the annihilation 
proceeds resonantly. 


This feature can be understood as follows. In Solution I, the heavy CP-odd Higgs 
boson is doublet-dominated with its mass usually at TeV scale. Then Eq.(2.3) in- 
dicates that the main contribution to the annihilation comes from the moderately 
light A,, which is singlet-dominated. With the formula presented in Eq.(A.20) and 
Us = u/À Z 450GeV shown in the lower left panel of Fig.3, one can get 


2 2 
2 2 2,2 ( Mo 5 2 
Chogo CA, bb = À K (=) < (5) K (3.1) 
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This inequation means that the couplings involved in the annihilation are highly sup- 
pressed so that the process must proceed resonantly to ensure (0,;0)9 ~ 107 78em3/s. 
Moreover, our results indicate that the width of A; is very small, rA, < 10~?MeV. 
So as ma, approaches 2m,o, the denominator in Eq.(2.3) tends to vanish and a small 
k in Eq.(3.1) is then suffice to predict the right rate of the annihilation for the GCE. 
This character is illustrated in the upper left panel of Fig.3. In fact, a small « is also 


favored to predict light ¥? and A1, which can be seen from Eq.(A.3) and Eq.(A.8). 


e The parameter u is upper bounded by about 300GeV, which is shown in the lower 
panels of Fig.3. 


This feature is actually required by the DM relic density [69]. Generally speaking, in 
order to predict the measured Qh?, the velocity weighted cross section (ov) should 
be around the canonical value 3 x 10726 cm/s at freezing out (see for example 
points in Table 2). Since 2m.o/ma, > 1 in Solution I, (ov) for the annihilation 
xix? o AT > bb at present day is usually larger than that at freezing out due to 
the thermal broadening [117]. Since the dwarf galaxy measurements have required 
(agvo S 2 x 10-76 cm3/s (see Fig.2), new contributions such as those mediated by a 
Z boson or a CP-even Higgs boson must intervene for the DM annihilation in early 
Universe, and a moderately small u can accelerate the annihilation [69]. 


e Solution I suffers from severe fine tuning problem. Explicitly speaking, beside the cor- 
relation m4, ~ 2m<0, there exits another strong correlation observed in our analysis, 
which is given by 


51 — 475k, for 30GeV < mgo < 40GeV or 0.024 < k < 0.045, 
myo /GeV œ~ 4 374+ 325K, for 50GeV < mgo < 62GeV or 0.038 < « < 0.07, 
49 + 175x, for 63GeV < mgo < 70GeV or 0.08 < «x < 0.12. 


These correlations make Solution I in the NMSSM quite unnatural to explain the 
GCE. 


e We checked that Br(A, > yy) < 5 x 1074 so that the y-ray spectral line generated 
by XL — At — yy is suppressed. 
< 60GeV for most cases in Solution I, the SM-like Higgs boson Hə may 


decay into ¥9 pair. We checked that Br(H2 > ¥?x?) < 18%, which is required by 
the Higgs data at the LHC. 


e Since x? 


3.2 Solution II - the HA; annihilation channel 


Solution II is quite similar to the interpretations presented in [49, 50, 83, 93, 99, 100], 
which utilize the process YX > ¢1¢2 > fıfıf2f2 ($1 and ¢2 denote scalar or pseudoscalar 
particles, and fı and fz are SM fermions) for the GCE. These interpretations, as were 
emphasized by the proposers, can easily escape the constraints from DM detection experi- 
ments and have been paid more and more attention recently. 

The features of Solution II are as follows: 
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e The singlet-dominated particles satisfy 60GeV S mgo < 115GeV, 10GeV S ma, S 
110GeV, 60GeV S my, S 120GeV and 6 < 0.2, and for most samples there exist 
following relations my, > myo > ma,. Given « ~ 0.1 which is required to predict 
the right size of the annihilation ¥?¥? + HA, for the GCE (see below), the particle 
spectrum limits parameters such as À, u and A, in certain regions (see the expressions 
of the tree level masses in Appendix), which are given in Table 1, and also shown in 
Fig.3. 


Note that u is below about 200GeV. In this case, the higgsino-dominated neutralinos 
X? may decay dominantly into ¥?A1 instead of into ¥?Z since the kinematics is 
forbidden. In this case, the LHC search for electroweakinos by trilepton +E/7"%S 
signal is less efficient in ruling out the light higgsinos?. Also note that the parameters 
à and p are related by u/GeV ~ 60 + 260 A for \ varying from 0.2 to 0.6 (see lower 
left panel of Fig.2), which means that vs = u/A > 360GeV. This ensures that the 
expansions for the masses and couplings in Appendix by the power of Av/j are good 
approximations. 


e The s-channel contributions to the annihilation rate (7, A; vo in Eq.(2.5) are usually 
much smaller than those from the t/u channel, and among the t/u channel contribu- 
tions, the one induced by the exchange of X} is far dominant. As for the contributions 
induced by the two higgsino-like neutralinos, each of them may be sizable, but since 
they cancel each other, the net higgsino contribution is not important. These char- 
acters can be understood by the following approximations (see Eq.(A.20)) 


Av 
C4, 3959 C mr, x099 = 2in? (1 + rk 


= P sin? 8, for Higgsino — like x? and mgo < 0, 
Cuat m9 ~) ity? oo Sad i 0 
Ii ei ria sin’ 6, for Higgsino — like x; and mgo > 0, 


and by the fact that x ~ 0.1 is enough to predict the ¥? contributed (oy,4,v)o0 at 
the order of 10~*°cem?/s (see equation (3.20) in [93]). 


e Since ma, < 60GeV for most cases in Solution II (see upper right panel of Fig.2), 
the SM-like Higgs boson Hə may decay into A ,Aj, with a sizeable fraction. Given 
that A; decays dominantly into bb, this will result in 4b signal for the SM-like Higgs 
boson. We checked that Br( Hz — A141) < 24%, where the upper bound comes from 
the constraints of the LHC Higgs data. 


e Since a good fit to the GCE requires that HA , is produced close to threshold, the 


annihilation ¥?¥? — HA, will produce spectral line or box-shaped spectrum in 


y-ray [83, 93]. We checked that Br(H, > yy) < 1 x 1073 for most samples and 


In doing [105], we once confronted with the situation quite similar to what we are facing now. Our 
detailed simulation at that time indicated that the trilepton constraint on SUSY is very weak. Moreover, 
in comparison with the case discussed in [118], we find that our case is more difficult to detect since the 
signal is smaller. 
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Figure 4. The GCE y? as a function of DM mass for Solution III where DM annihilates mainly 
into WW final state. Here Ry, 4, denotes the branching ratio of the annihilation into H1 A, final 
state. 


Br(A, > yy) < 4x 1074 for all samples. So current results of the Fermi-LAT search 
for spectral lines [119] can not impose tight limit on Solution II (see [83] for a detailed 
discussion). 


3.3 Solution III - the WtW- annihilation channel 


In general, the pure annihilation ¥?%? — W*W is unable to explain the GCE quite well 
(73, 79], but if it mixes sizably with other annihilation channels, the generated spectrum 
may be improved significantly to account for the GCE. Solution II in the NMSSM belongs 
to this case. 

Solution II has the following features: 


e The W pair must be produced close to threshold to account for the GCE, which 
means 85GeV < mgo < 100GeV (see right panels of Fig.3). 


e From the expression of (owwv)o in Eq. (2.4), one can learn that, if the wino is de- 
coupled, the annihilation rate is determined by the higgsino-dominated chargino. In 
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this case, we have 


fin S-AN 14 = g 8 = VPasin EE, 


fir > -3M ~ -V2gsin p. (3.2) 
In getting these expressions, we note v, = u/A 2 400GeV (see lower left panel of 
Fig.2), and expand Nı3 and Nj, in terms of Av/ (see Appendix). We also use the 
approximation myo ~ 2ku/à. Then (owwv)o ~ 1077%em3/s and ms xo ~ 90GeV 
limit tightly the es of the parameters A, k and u, which are iao in Table 1 
and Fig.3. 


Note in Solution II, the parameter u, or equivalently the masses for the higgsino- 
dominated chargino and neutralinos, is less than about 150GeV. Since the splitting 
between u and mgo is less than about 50GeV, such a low value of u is still allowed 
by the LHC search for SUSY (see footnote 2 in our discussion on Solution I). 


e The upper left panel of Fig.3 indicates that the parameters À and « are correlated by 
k œ 0.03 + 0.3å, for 0.2 < A< 0.4. (3.3) 


As a result, we have mgo ~ 2u/3. 


e As we emphasized before, the annihilation X9} — WtW~ must mix sizably with the 
annihilation ¥?X? — HA to explain the GCE. This, in return, requires appropriate 
masses for Hı and A; to improve the y-ray spectrum generated by the WW state. 
In Fig.4, we plot the GCE x? as a function of DM mass in Solution III with different 
colors denoting the branching ratio of the DM annihilation into H,A,. This figure 
indicates that, with the increase of the branching ratio, the GCE y? tends to decrease. 


4 GCE solutions with H, being the SM-like Higgs boson 


In this section, we investigate the GCE solutions for the case that DM is singlino-dominated, 
and meanwhile H; acts as the SM-like Higgs boson. We carry out our study in a way similar 
to what we did in Section III. 


In Fig.5 we project the solutions on (av)9 — mo plane (upper panel) and x&cop — M0 


plane (lower panel). For solutions marked by green lozenge, DM annihilates with the ie 
branching ratio into Hı A1, while for those marked by red pentastar, DM annihilates mainly 
into H2Aı. In the following, we call these two kinds of solutions Solution IV and Solution V 
respectively. Fig.5 then indicates that, for 70GeV < mgo < 87GeV, Solution IV can explain 
the GCE quite well with the best explanation coming from mgo ~ 81GeV with Noe gp œ 21.4 
(corresponding to a p-value of 0.55), and for 80GeV < mo < 130GeV, Solution V is good 
in accounting for the GCE with the best explanation locating at mgo ~ 85GeV with 
XŽcr % 21.6 and a p-value of 0.54. 

Compared with the case that Hə acts as the SM-like Higgs boson, we find that it is 


more difficult to get the GCE solutions if Hı corresponds to the SM-like Higgs boson. One 
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Figure 5. Similar to Fig.2, but showing the solutions for the case that Hı acts as the SM-like 
Higgs boson. For these solutions, DM may mainly annihilate into Hı A; final state (called Solution 
IV in our discussion) or into H2A, final state (Solution V). 


important reason is that the spectrum of the singlet-dominated particles for Solution IV 
and V has non-trivial requirements on the NMSSM parameters, which can not be easily 
satisfied due to the structure of the NMSSM itself. A good example about this argument is 
that we do not find any solutions where DM mainly annihilates into bb. This is due to the 
fact that, given a singlino-dominated DM with 30GeV < myo < 70GeV and meanwhile a 
singlet-dominated A, satisfying m4, ~ 2mzo, the singlet-dominated CP-even Higgs boson 
is usually lighter than the SM-like Higgs boson [69]. 

In Table 2, we present detailed information for benchmark points P4, P5, P6 and 
P7 with points P4 and P5 belonging to Solution IV and points P6 and P7 belonging to 
Solution V. This table shows that for 80GeV < m<0 < 86GeV, DM may annihilate into 
HA, and H24; states with comparable rates to explain the GCE (see points P5 and P6), 
while for mgo ~ 78GeV (myo ~ 120GeV), the sole annihilation channel ¥9¥? > Hi Ai 
(XLX — H2A1) can be responsible for the GCE, see point P4 (P7). 

In Fig.6, we show the correlations among the parameters A, K, H, mgo and ma,. 
This figure is supplement to Table 1, and as we will show below, it is helpful for our 
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Figure 6. Similar to Fig.3, but showing the correlations for Solution IV and Solution V. 


understanding on Solution IV and V. 


4.1 Solution IV - the A, A, annihilation channel 


Solution IV has the following features: 


e The HA, state must be produced close to threshold to explain the GCE, which is 


reflected by 6 < 0.1 from our results. 


e The favored spectrum for the singlet-dominated particles is 71GeV < mgo < 87GeV, 


10GeV S ma, < 40GeV and 126GeV S mp, 


N N 


< 142GeV. Given & ~ 0.12 which is 


required to predict the right size of the annihilation %9} — HA, for the GCE (see 
below), this spectrum limits parameters such as A, «, 4 and A, within rather narrow 


ranges, which are given in Table 1 and also shown in Fig.6. 


Compared with Solution I, we find in Solution IV that, in order to predict a heavier 
singlet-dominated CP-even Higgs boson, the parameter u usually takes a larger value, 
220GeV S uw S 270GeV. As a result, A must exceed about 0.6, which can be inferred 
from the relation mgo ~ 2ku/A ~ 2 x 0.12 x u/A ~ 80GeV. This relation also 


= 17— 


suggests that vs = u/A 2 300GeV or Av/u < 0.6, which makes the expansions listed 
in Appendix feasible. 


e Similar to Solution II, the s-channel contributions to the annihilation rate (oH, Avo 
in Eq.(2.5) are usually much smaller than those from the t/u channel, and among 
the t/u channel contributions, the one induced by the exchange of X} is dominant. 
However, since Hı now is the SM-like Higgs boson (instead of a singlet-dominated 
particle in Solution II), there still exists a slight difference between the two solutions. 
Explicitly speaking, we find that each higgsino contribution to the annihilation is 
comparable in magnitude with the ¥? contribution, but since the two higgsino con- 
tributions cancel each other, the total higgsino contribution is small. These features 
can be explained by the following formula (see Eq.(A.20)) 


K 


sin? 8, (4.1) 


Av 
“2 SD pp Be eB 
C41 5999 CH, 99 Aik 7 sinf B ~ —8ik 


mie 


c c mwa sin? 8, for Higgsino — like X and meso < 0, 
A131% MRI ~ -iya sin? 8, for Higgsino-— like X? and mgo > 0, 
and also by comparing Eq.(4.1) with equation (3.20) in [93] to conclude that « ~ 0.12 
is enough to predict the ¥? contributed (cy, 4,V)o at the order of 107? cm’ /s. 


e Since m4, < 40GeV for all cases in Solution IV (see upper right panel of Fig.6), the 


SM-like Higgs boson Hı will decay into A;A;. We checked that Br(H, > A,A1) S 
24% as required by the Higgs data. 


e We also checked that A; — bb is the dominant decay mode of A1, and Br(A; > 
yy) <5 x 1075 for all samples. 


4.2 Solution V - the HA, annihilation channel 


Since Hə in Solution V is singlet dominated, the features of Solution V should be similar 
to those of Solution II. The differences mainly come from the following aspects: 


e The spectrum of the singlet dominated particles. In Solution V, the favored spectrum 
is 80GeV < mgo S 130GeV, 18GeV < ma, £ 100GeV and 125GeV S mp, S 
146GeV with ma, < mgo < MH, and 6 < 0.1. Corresponding to such a spectrum, 
the parameter space of Solution V differs greatly from that of Solution II, which can 


be seen from Table 1 and also from Fig.6. 


e The phenomenology of some relevant particles. For example, in both Solution IV and 
Solution V, the favored value of u is uplifted in comparison with that in Solution II. 
As a result, the higgsino-dominated neutralinos may decay into Z%9, which makes 
them to be potentially detected at 14-TeV LHC by trilepton +E7"*S signals [69]. 
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Figure 7. Spin-independent (SI) and Spin-dependent (SD) cross sections for DM-nucleon scattering 
as a function of DM mass. Solutions in this figure are taken from Fig.2 and Fig.5 with the same 
symbolic conventions. 


5 Explore the GCE solutions in future DM experiments 


In this section we investigate to what extent the GCE solutions will be explored in future 
DM direct detection experiments such as XENON-1T and LUX experiments [120], which 
will improve current experimental sensitivities to DM-nucleon scattering cross sections by 


SI SD 


up to three orders. In Fig.7, we project our solutions on m o?’ and mgo — o 


ppe 
planes with a and rte denoting the spin-independent (SI) ad TOO ca (SD) 
cross sections respectively. The left panels in the figure are the results for the case that 
Hə acts as the SM-like Higgs boson, and the right panels are those for the case that Hı 
corresponds to the SM-like Higgs boson. The dotted lines, solid lines, dashed lines and 
dash dotted lines are the sensitivities to the cross sections set by the XENON-100, LUX, 
XENON-1T and LZ experiments respectively. Note that so far the XENON-100 experiment 
has imposed constraints on both SI and SD cross sections, while the LUX experiment only 
obtained limits on the SI cross section. 

For a in the Hə case, we can see from Fig.7 that the future XENON-1T experiment 
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is able to probe a large portion of the GCE solutions, and the LZ experiment can test even 
more solutions. Anyhow, there still exist some solutions remaining untouched by these 
future experiments. This conclusion can be understood as follows. In the NMSSM after 


considering the current experimental constraints on sfermion masses, the main contribution 


I 


to os comes from the t-channel process mediated by the CP-even Higgses H1,2. In this 


case, the Wilson coefficient fy, for the operator X9 X}g;q; is given by [121] 


Pa Ch909 C Hiqi C H2399 C Hoqiqi 
qi — I 


(5.1) 


2 2 
2M 2M 


where Cy, gog ~ —v2k(1 + 2\v/p) and Cm S 2V2kAv/u, which are given by 
Eq.(A.20) and Eq.(A.21) respectively. Then Eq. (5.1) indicates that, if x is small or if 
there exists a strong cancelation between the two terms, fg, or equivalently the SI cross 
section will be suppressed. We numerically checked that the untouched solutions has either 
of the two characteristics. 

On the other hand, the story for oo” in the Hə case is quite different. From the 
lower left panel of Fig.7 we can see that the future XENON-1T experiment can test almost 
all of the GCE solutions, let alone the more sensitive LZ experiment. The underlying 
reason is that in the NMSSM with heavy sfermions, the SD cross section gets contribution 
mainly from the ¢-channel Z-mediated diagram. As a result, the size of the cross section 
is determined by the Z7¥?¥? coupling, which is given by 


22 AK? mz Hv 
“epee Nae Wee ~ 1 KLJ; 5.2 
J9Zzz9%9 Jon 13 i4) Vw ye ( y7) Vw ye ( ye ) ( ) 


In getting this expression, we have used the approximations for N13, Ni4 and mgo. Then 
from the results presented in Fig.3, one can infer that except for some rare cases of Solution 
I, the SD cross section is not suppressed too much. 

In a similar way, one can analyze the results for the Hı case, which are shown on the 
right panels of Fig.7. For example, the upper right panel indicates that the SI cross sections 
in Solution IV and Solution V are usually larger than 10~!9pb. This may be understood by 
a weak cancelation between the two terms in Eq.(5.1). Compared with the Hə case, both 
the SI cross section and the SD cross section in the H; case are large and consequently, all 
the solutions will be tested by XENON-1T experiment. 

In principle, the GCE solutions in the NMSSM may also be tested by electroweakino 
production processes at the LHC [69]. We will discuss such an issue in our forthcoming 
work. 


6 Summary 


In this work, we took into account the recently reported foreground and background un- 
certainties for the GCE and investigated its explanation by DM annihilation in the frame- 
work of the NMSSM. We carried out a sophisticated scan over the NMSSM parameter 
space by considering various experimental constraints such as the Higgs data, B—physics 
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observables, DM relic density, LUX experiment and the dSphs constraints. Then for each 
surviving parameter point we performed a fit to the GCE spectrum by using the corre- 
lation matrix that incorporated both the statistical and systematic uncertainties of the 
measured excess. Our results indicate that due to the introduction of the gauge singlet 
Higgs superfield, the NMSSM with a singlino-dominated DM has multiple DM annihilation 
channels that are able to explain the GCE quite well, and all of these explanations require 
the singlet-dominated particles (including one neutralino, one CP-even and one CP-odd 
Higgs bosons) to be moderately light. We also discussed to what extent the future DM 
direct detection experiments can explore the GCE solutions, and we conclude that the 
XENON-1T experiment is very promising in testing nearly all the solutions. 

When choosing the scenario of particle spectrum, we focused on a singlino-dominated 
DM and considered the cases that either Hə or Hı acts as the SM-like Higgs boson. For 
the popular situation that Hə corresponds to the SM-like Higgs, we have the following 
observations on the GCE solutions: 


e The pure DM annihilation channel ¥?x? — bb or X9X? —> Hi Aı can provide a good 
fit to the GCE spectrum, while the channel ¥?x? ~ WWT must mix sizeably with 
the channel %9? > Hı Aı to account for the GCE. 


e For the annihilation ¥?X%? —> bb, DM mass is now allowed in the range from 30GeV 
to 70GeV which is much wider than before. With the help of an appropriate s- 
channel resonance, the singlet trilinear self-coupling parameter «x can be as low as 
0.02 to explain the GCE. Moreover, the higgsino mass parameter u is upper bounded 
by about 300GeV to ensure a correct DM relic density. Since there exist strong 
correlations between independent parameters, such an explanation suffers from a fine 


tuning problem, which is usually less than 1%. 


e The annihilation ¥?X%? > HA may provide a better explanation than the channel 
Xx? — bb when Hj A, is produced close to threshold, and the best interpretation 
corresponds to a p-value of 0.55. In this kind of explanation, the singlet-dominated 
particles must satisfy 60GeV < mgo < 115GeV, 10GeV S ma, < 110GeV, 60GeV < 
my, S 120GeV and 6 < 0.2. This imposes non-trivial constraints on the NMSSM 
parameters, especially that u must be less than about 200 GeV. Among various 
contributions to the annihilation, the dominant one comes from the ¥-contributed 
t/u channel diagrams, in which the parameter « plays an important role in deciding 


the annihilation rate. 


e Apart from the necessary mixing with the H, A final states, WtW- pair in the 
annihilation ¥?x? — WHW- must be produced close to threshold to account for the 
GCE. A small p less than about 150 GeV is necessary to increase the annihilation rate 
through the t/u-channel contributions induced by a higgsino-dominated chargino. 
The LHC search for trilepton + Eipiss signal can not exclude such a possibility since 
the electroweakino production rates at the LHC are relatively low, and meanwhile 
since the splitting between u and mgo is compressed. 
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e The detection of spin-independent scattering in the future XENON-1T and LUX 
experiments are able to cover a large portion of the GCE-favored parameter space, 
while the spin-dependent detection have a stronger potential to test nearly all of the 
relevant parameter region. 


As for the case that Hı acts as the SM-like Higgs boson, the features of the GCE 
solutions are quite different, which are as follows: 


e In comparison with the Hə case, it is difficult to find GCE solutions when H; corre- 
sponds to the SM-like Higgs boson, and especially we did not find any solution that 
DM annihilates mainly into bb. The reason is, assuming Hı to be the SM-like Higgs 
boson, there must exist sizeable mass splittings among the light singlet-dominated 
particles to explain the GCE, which is difficult to realize in the NMSSM due to the 
theoretical structure itself. 


e For 80GeV < Mo < 86GeV, DM may annihilate into H,A, and H2A, states with 


1 
comparable rates to explain the GCE, while for mgo S 80GeV (myo = 100GeV), 
the sole annihilation ¥9¥9 > Hi Ai (X?X? —> H2A1) can be responsible for the GCE. 
For all these solutions, the singlet-dominated particle Hə and the parameter u must 


satisfy 125GeV S my, S< 145GeV and 210GeV S ps S< 270GeV. 


e Both the spin-independent and spin-dependent detection in the future XENON-1T 
experiment have a great potential to test the relevant parameter space. 


Before we end our discussion, we would like to comment briefly on the interpretation 
of the GCE with a bino-like DM. Like the singlino-dominated DM case, a light A, with 
mass below about 140GeV is necessary for such a work, and this A; prefers to be singlet- 
dominated °. The difference is that, for the bino-like DM case, the interaction of the DM 
with A, is relatively small and consequently the annihilation ¥?¥? > H;Aj can not explain 
the GCE any more due to its rather low annihilation rate. Also due to the suppressed 
interaction, ma, must be closer to 2M0 for the annihilation X9? — bb to account for the 
GCE, and thus the theory has to be tuned in a more elaborated way. Our sophisticated 


scan over the relevant NMSSM parameter space verified these conclusions. 


3In the NMSSM, a light Ai with mass below about 100GeV may have a large doublet component if the 
elements of the CP-odd Higgs mass matrix satisfy M99 > M19 > Moi (see benchmark points P3 
and P4 in [122]) or M? a22 ~ Moi ~ M? 12 (see the point presented in Table 2 of [123]). In either case, 
ma, should be significantly smaller than Mmy+ to escape experimental constraints. Previous studies have 


suggested that a light doublet-dominated A; might also explain the galactic center excess. However, due 
to the requirements on the elements this scenario occurs only in specific portions of the parameter space 
and is significantly more experimentally constrained than those we considered. In fact, in our scans for the 
GCE we did not find any parameter points with the doublet component of the light A; exceeding 0.1. In 
summary, a light doublet-dominated A; may exist, as suggested by e.g. Ref. [123], but it is fair to say that 
without a very delicate parameter tuning, it is difficult to obtain in explaining the GCE, especially when 
one considers more constraints than previous literatures. About this conclusion, we thank the authors of 
[123] for helpful discussion. 
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A Properties of the singlet-dominated particles 


In this appendix, we present some analytic expressions for the masses and couplings of 
the singlet dominated particles, such as ¥? and A, in the NMSSM. These expressions are 
obtained by diagonalizing the mass matrices of the particles (like done in [57]), and are good 
approximations in certain cases. They are helpful in understanding the results presented 
in this work. In the following, we will follow notations and conventions consistent with 


[104] for the Z3 NMSSM. 


A.1 Neutralino masses and mixings 


In the basis Y° = (—ià1, —iA3, Y9, Y}, ws), the neutralino mass matrix is: 


v2 V2 
Mo oe = 0 
M= 0 =u —Avy |: (A.1) 
0 —Avg 
pau 


If the bino and wino fields are decoupled, the mass eigenstates of the neutralinos can be 
approximated by 


E & Migb9 + Nia? + Mises, 
Xo = Nigh + Nayt + Nisvs, (A.2) 


where ¥? denotes the lightest neutralino with Yg field as its dominant component in this 
work, and ¥? represents a higgsino-like neutralino. 

In the limit of |u| > Av, 1 > k/À and tan > 1, one can expand the neutralino 
masses and N;; by the power of Av/u = v/vs to get the following approximations: 


mo ~ Sut M (usin 28 = a) 

a a = ae cos 8 (tan Braye — n) x — sin, 

ia ua aa 

Ms = (1+4 a L Na ~ ZSan(inga)Almse), 

Na = —ZeSan(u)O(mse), Nam- E sgn). (A3) 


— 23 — 


In above expressions, the Sgn and 0 functions are defined by 


smod yee, Wl, deen AM 


Likewise, one may consider the case that the wino and the singlino fields decouple. In 
this case, the mass eigenstates of the neutralinos is approximated by 
~O 0 0 
Xi S Nul idi) + Nispa + Nady, 
XP ~ Na(—iA1) + Niph + Nid. (A.5) 


In the limit of tan 8 > 1, |u| > gavy and |u| > Mı, we have the following approxima- 
tions: 


De 24 

mgo % Mı — ae aan 6 + Mı), 

1 p 
Ni3 _ mzsin ðw ing 
Nu 
N: ind M 

He A cos 3(1 + tan B—), 
Nu H u 

Nia N 1 
Nu = (1 ae N2 + N2 x1, Ni% Pe rm g9 
1 1 mz sinw sin 8 

Nabio == Dane): DLS 7 Sgn(u)O(mgo). (A-6) 


A.2 CP-odd Higgs mass matrix 


In the (A, Sr) “interaction” basis, the mass matrix for CP-odd Higgs bosons is given by 


ma . 3 
m, Avla sin 26 — 4) 


ee mi 3K) _ 3KA 
2v2 sin 26 (72 sin 26 + 3%) — 34st 


(A.7) 


In the case of ma > max(v,|Ax|, lul), K/À < 1 and tan 8 > 1, the lighter CP-odd scalar 
A, is singlet dominated with its squared mass given by 


3K A, LL 
A 


This approximation indicates that, without considering the radiative corrections, the singlet- 


9 
Må 7 au sin 28 — 5 (A.8) 


dominated CP-odd scalar mass is determined by the parameters A, x, u as well as A,. The 
components of A, can be written as 


P A 2 d 
ee & sin 28 St) ~ ——" sin 26, 
PAy,sr må \ 2 2u 


p? -1/2 
Ai,A 
PAS = (: BE a=) x1, (A.9) 


? 


where P4, ,a is the active component and P,4,,5, is the singlet component of the Aj. 


— 24 — 


A.3 CP-even Higgs mass matrix 


In the basis (S1 = cos Bhu — sin Bha, S2 = sin Bhu + cos Bha, S3 = hg), the mass matrix el- 
ements for the CP-even scalars are [107] 


Mi, = Må + (mz — d*v”) sin? 28, (A.10) 
1 
Mig = —5(mz — Xv’) sin 48, (A.11) 
2p, À 
M?, = —(M3 sin 28 + a K cos 2, (A.12) 
Mi, = m3 cos? 28 + X°? sin? 28, (A.13) 
Ma, sin 26 Kj, 
23 = : A14 
M53 = 2Auv |1— ( Dy ) zy sin 26|, ( ) 
May = e | “(An + wh) = pase? sin 28, (A.15) 


where Sz is nothing but the Higgs field in the SM, M3, is its mass at tree level without 
considering the mixing among Sj, and the second term \?v? sin? 23 in M3, originates from 
the coupling AH, - Aa in the superpotential. 

The mass eigenstates H; are defined by 


H; = Va S1 + Vi2Se + Vis S3, (A.16) 


where V is the rotation matrix to diagonalize the mass matrix. For the S2-dominated mass 
eigenstate H;, current Higgs data have required it to be highly SM-like, i.e. Vj1,Vj3< 1, 
so in the case of the hierarchy structure M7, > maz( M3s, M33), M33 K |M3. — M33]. 
If we decouple the MSSM-like heavy Higgs, S1, from the other two, the 2 x 2 reduced mass 
matrix in the (S2, 83) basis is given by [107]: 


Me Mas i A.17 
SaR = M3, — “mA tia sin? 48 — © vals V? cos? 28 — Amv? cos? 28 sin 28 ) 


The (2,2) element of the reduced (2 x 2) matrix, which in the limit of zero-mixing with the 
other Higgs should give singlet scalar mass in the Z3 NMSSM, is given by: 


4 Pm 
M3,s;(2, 2) = z (4. + re) poo (1 — cos? 28) sin? 28 


A À 4u? 
2,2,7,2 
1 
-EPP cos? 28 — — ene (2 cos? 28 + 1) sin 2. (A.18) 
må 2 
Setting M3, ~ 0, i.e. m4, = (1 — 3x sin 26), and taking tan 6 > 1, we have 
K Ak 
Mi, ¢, (2,2) © T (4, £ =) (A.19) 


This approximation indicates again that, without considering the radiative corrections, the 
singlet-dominated CP-even scalar mass is determined by the parameters \, <, y and A,. 
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A.4 Some properties of the singlet-dominated particles 


With the assumptions that M4 > maz(|ul,|Ax|), |u| > Av, tan’ > 1 and K/A < 1, 
one can approximate the masses and couplings of the singlet dominated particles, such 
as ¥? and Aj, by simple analytic expressions [57]. In the following, we list some of the 
coupling expressions used in our discussion, which are denoted by Cxyz hereafter. These 
expressions are actually expand the corresponding exact ones by the power of Av/,1. 


Cre NON p, dee imp Au 
! vw V2v u 
—iV/2n(1 +242) for y= 1; 
Car x09 ~ i — Ae Sgn(mzo)O(mzo) for higgsino — like X, 
J I 
CHa = iC A, if H; is singlet dominated. (A.20) 


22K Ae forj=1, 
1Xj 3 Sgn(mzo )O(m<o) for higgsino — like X;, 
J I 
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